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Abstract 

The RHIC experiments have measured the nuclear modification factor Raa of 
non-photonic electrons in Au+Au collisions at ^Jsnn = 200 GeV. This Raa 
exhibits a large suppression for pt > 2 GeV/c which is commonly attributed 
to heavy-quark energy loss. It is expected that the heavy-quark radiative 
energy loss is smaller than the light quark one because of the so-called dead- 
cone effect. An enhancement of the charm baryon yield with respect to the 
charm meson yield, as it is observed for light and strange hadrons, can explain 
part of the suppression. This phenomenon has been put forward in a previous 
work. We present in this paper a more complete study based on a detailed 
simulation which includes electrons from charm and bottom decay, charm and 
bottom quark realistic energy loss as well as a more realistic modeling of the 
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h-c/D enhancement. We show that a h-c/D ratio close to unity, as observed 
for hght and strange quarks, could explain 20 — 25% of the suppression of non- 
photonic electrons in central Au+Au collisions. This effect remains significant 
at relatively high non-photonic electron transverse momenta of 8 — 9 GeV/c. 
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One of the most robust experimental evidence for the creation of a new state of matter in 
heavy ion coUisions at the Relativistic Heavy Ion CoUider (RHIC) is the large suppression of 
light hadrons at high transverse momentum (pt) [!]• This phenomenon is well reproduced by 
models which take into account the radiative energy loss of high pt light quarks and gluons 
propagating through a dense medium of colored quarks and gluons [2] . Further insights into 
the underlying mechanism can be obtained from the study of heavy hadrons. In contrast 
to intermediate-pt light hadrons which are predominantly produced by gluon fragmentation, 
charm and bottom hadrons originate from the fragmentation of heavy quarks. Quarks 
are supposed to lose less energy than gluons in the medium due to a smaller color charge 
coupling. In addition, radiative energy loss was predicted to be smaller for heavy quarks 
as compared to light quarks because of the so-called "dead-cone" effect which limits the 
medium induced radiative energy loss at forward angles [3]. Surprisingly, recent data from 
the PHENIX and the STAR collaborations in Au-^Au coUisions at y^Sivjv = 200 GeV show 
that the quenching of heavy quarks, as studied indirectly via the so-called non-photonic 
electrons^, is stronger than theoretical expectations [4-6] and is as large as that of light 
mesons. Reconciling these data with model predictions is a real challenge which triggers 
a lot of theoretical activities nowadays. Only models which assume a very large medium 
opacity [7], an additional coUisional energy loss [8] or effective energy loss from multiple 



^In contrast to light hadrons, the heavy flavor quenching is, so far, not measured experimentally 
through identified hadrons, but in an inclusive way via the nuclear modification factor {Raa) 
of non-photonic electrons. The latter is obtained from the pt distributions of electrons (after 
subtraction of Dalitz-decay electrons from light hadrons and photon-conversion electrons) in AA 
collisions {dN^^/dpt) and in pp collisions (dN^/dpt) as: 

dN%Jdpt 

< N^J > dNe^/dp, 

where < > is the average number of nucleon-nucleon collisions corresponding to a given 

centrality class. 



fragmentations and dissociations of heavy quarks and mesons {D and B) in the medium [9] 
can describe, with a relatively good agreement, the data (for a recent review, see [10]). 

In this paper we investigate the possibility that part of the strong suppression of non- 
photonic electrons might be due to another source of electrons, namely charmed baryons. 
Indeed, whereas hght mesons are largely suppressed in heavy ion coUisions at RHIC, the 
suppression of non-strange and strange baryons is observed to be much less in the inter- 
mediate pt range (2 < pt < 4 GeV/c) [11]. This is commonly referred to as the anoma- 
lous baryon/meson enhancement. This anomalous baryon/meson enhancement is relatively 
well understood in the framework of the recombination model which assumes that, at low 
and intermediate pt, hadronization occurs via the coalescence of "free" quarks (and anti- 
quarks) [12]. An anomalous baryon/meson enhancement for charm hadrons leads naturally 
to a non-photonic electron Raa smaller than one. This is mostly due to a smaller semi- 
leptonic decay branching ratio of charm baryons (Ac) as compared to charm mesons (see 
Tab. I). As a consequence, part of the experimentally measured Raa of non- photonic elec- 
trons should not be attributed to energy loss. We show that the AJ D enhancement can 
explain up to 25% of the non-photonic electron suppression data measured by the PHENIX 
collaboration in Au-|-Au coUisions at i/sjvjv = 200 GeV [13]. 

The main assumption we put forward is that, in a deconfined medium, charm baryon 
production is enhanced relative to charm meson production, as compared to the vacuum. 
This assumption is qualitatively justified in the framework of the recombination model. 
Although this model docs not provide detailed predictions on charm hadron production yet, 
it successfully describes the (non-charm) baryon/meson enhancement measured in Au-|-Au 
collisions at ^/snn — 200 GeV. A relatively good agreement is obtained not only for the 
light hadron ratio p/ir^, but also for heavier hadron ratios such as A/K^ and ^2/0 [14]. 
Extrapolating these results to charm hadrons is not straightforward because the mass of 
the charm quark is much larger than that of light and strange quarks. The consequences 
are threefold as far as the recombination mechanism is concerned. First, whereas the pt 
of a light baryon (meson) amounts to 3 (2) times the initial pt of its valence quarks, the 
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Pt of a (single) charm baryon or a charm meson is hkely to be very close to that of the 
charm quark. Secondly, considering a light quark and a heavy quark with the same velocity 
(which is the essential requirement for the coalescence process to take place [15]), the heavy 
quark momentum is much larger than that of light partons. As a consequence, one can 
expect the enhancement of the charm baryon/meson ratio to appear at higher p^ than 
that of non-charm hadrons. The recombination model indeed predicts, for non-strange 
and strange hadrons, that the heavier the hadron, the larger the pt of the baryon/meson 
enhancement. This has been observed for non-strange and strange baryon/meson ratios by 
the STAR collaboration [14]. Finally, the fragmentation time of heavy quarks is small as 
compared to hght quarks. According to [9], the formation time of a 10 GeV/c pion, D 
meson and B meson is 20, 1.5 and 0.4 fm/c, respectively and it is as small as ~ 3 fm/c 
for a Ac with p^ = 20 — 30 GeV/c. Due to these considerations, it is obvious that the 
baryon/meson enhancement for non-charm hadrons and charm hadrons can be significantly 
different. In the following, we only assume that, in view of experimental results on the 
baryon/meson enhancement for non-strange and strange hadrons, a similar enhancement 
is a priori conceivable for charm hadrons^. Remarkably, such an enhancement has strong 
implications on the nuclear modification factor of non-photonic electrons. It leads to a 
decrease of the yield of non-photonic electrons in A + A collisions because, as shown in 
Tab. I, the inclusive semi-leptonic decay branching ratio of charm baryons is smaller than 
that of charm mesons. Therefore, the nuclear modification factor of non-photonic electrons 
should decrease as well. This can be easily illustrated in the following way. Assuming that 
charm production scales with the number of binary collisions (i.e. Raa = 1 in absence 
of medium effects) and that the relative yields of D mesons are the same in pp and in 
A + A collisions, a pt integrated Raa can be calculated for different C enhancement factors. 



very recent theoretical study in the framework of the recombination model confirms this 
assumption [16]. 
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with 



according to 



N - IN ^ + (l\ 

Ac,Ac/ D ^ ^ ^ ^ ^ K ) 



where 



M /AT ^ (^Ac,Ac/ ^d)BRa,,Ac /ON 

e-Aj e^D ^Nd± / N D)BRn± + {N / N d)BRj,o + {N ^± / N d)BR^± ^> 
N is the charm hadron yield and BR is the hadron semi-leptonic decay branching ra- 
tio. According to Tab. I, N^^^j^JNd = 7.3%, Nd±/Nd = 21%, Nj^o^d^/Nd = 67%, and 
Nj-)±/Nd = 12% such that Ng^\^/Ng^D = 3.63%. Therefore, an enhancement factor C of 
12 leads to a non-photonic electron Raa of 0.79 ±0.07 and in the extreme case of an infinite 
enhancement, the non-photonic electron Raa reaches 0.51. 

The above idea has already been proposed in [17]. Before going to our simulation results, 
we present in Tab. II the main differences between our approach and the one of ref. [17]. The 
choice of a Gaussian shape for the pt dependence of the Ac/ D ratio in Au-|-Au collisions is 
motivated by results from the coalescence model for heavy quarks [19]. For pp collisions, we 
use the predictions from PYTHIA and not the shape from the measured A/K^ ratio since 
experimental results from the STAR collaboration indicate a strong mass dependence of 
baryon/mcson ratios [14]. These assumptions lead to a significant difference in the maximum 
the Ac/D ratio and in its location in pi, as it is reported later. Finally, the present work 
includes a more realistic treatment of the heavy-quark energy loss as well as the contribution 
of electrons from bottom decay which is ignored in [17]. 

Our simulation framework is based on the PYTHIA-6.152 event generator [20]. The 
PYTHIA input parameters were first tuned according to [21] and the PHENIX acceptance 
cut {\ri\ < 0.35) was applied in order to correctly reproduce the pt distribution of non- 
photonic electrons measured in pp collisions at ^/s — 200 GeV [13]. As it can be seen in 
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Fig. 1, the agreement between the simulation and the data is rather good except in the 
high pt region where the simulation under-predicts the data. The result of the simulation is 
also compared to FONLL (Fixed Order Next to Leading Log) predictions [22]. As already 
observed in [13], the PHENIX data is in agreement with FONLL within the theoretical 
uncertainties. 

Table I shows that the Ac/ D ratio amounts to 7.3% (in 47r) which translates to 3.63% 
after convolution of the species yields with their corresponding semi-leptonic decay branching 
ratio. On the other hand, Fig. 2 shows that, in the pt > 2 GeV/c region of interest discussed 
hereafter, this ratio is even smaller (~ 1.5%) because the decay electron spectrum of Ac is 
softer than that of D mesons. This leads to an additional suppression of the non-photonic 
electron yield at intermediate Pf 

As stated above, the non-photonic electron pt distribution in Au-|-Au collisions has been 
evaluated after considering a enhancement whose shape is, according to the predictions 
of the coalescence model [19], assumed to be a Gaussian versus Pf It has the following 
parameters. Mean: 5 GeV/c, constant: ~ 0.9 and sigma: 2.9 GeV/c. The constant of 
0.9 is obtained from N^^-^JNd x C with Nf^^-^jNo = 7.3% (Tab. I) and C = 12. Such 
an enhancement factor C = 12 is justified since the resulting h-c/D ratio of ~ 0.9 is of the 
same order of magnitude as the non-strange and strange baryon/meson ratios measured by 
the STAR collaboration [14]. In contrast, the corresponding (enhanced) h-c/D ratio is in 
ref. [17] located at lower pt and its maximum is close to 1.7. The enhancement is applied such 
that the pt-differential charm cross-section is conserved. The latter is an arbitrary choice 
that could be justified since most of the charm hadron transverse momentum is given by the 
charm quark whatever, baryon or meson, this hadron is. We finally compute the Raa ratio 
from the non-photonic electron pt spectra assuming that the only medium induced effect is 
the Ac/ D enhancement. The results are shown in Fig. 3 together with the PHENIX data. 
Note that at this step only electrons from charm decay are considered and heavy-quark 
energy loss is neglected. The simulated Raa ratio is shown only for pt > 2 GeV/c since 
shadowing might play a role and has not been considered in the simulation. One can see 
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that the Ac/-D ratio close to unity in central coUisions at pt = 5 GeV/c can already explain 
~ 40% of the suppression of non-photonic electrons in the 2 — 4 GeV/c pt range. Even in 
the high pt region (8 — 9 GeV/c) the Ac/ D enhancement results in a significant suppression 
of non-photonic electrons. 

In the next step c quark radiative and coUisional energy loss is included in the simulation. 
This is achieved by a convolution of our non-photonic electron pt spectra with the differential 
suppression factors taken from [23]. It is shown in Fig. 4 that the relative suppression 
originating from the Ac/ D enhancement is about the same amplitude than the one from the 
charm coUisional energy loss and represents about 36(20)% of the observed suppression at 
Pt — 4(9) GeV/c. In contrast, the suppression reported in [17] is less than 20% in the pt 
range 2 — 5 GeV/c and becomes negligible for Pt > 5 GeV/c. 

Finally the bottom contribution is added in the simulation. This obviously reduces the 
suppression of the sum of non-photonic electrons because b quarks are supposed to lose 
less energy than c quarks. However, the relative contribution of c and b quarks to the 
total non-photonic electron yield is not well known. According to FONLL predictions, the 
crossing point (p^^) between charm and bottom electron decay pt spectra is expected to be 
located in the range 2.5 < pt < 10.5 GeV/c. Therefore we have considered two scenarios to 
include the bottom contribution: a crossing point in the central value predicted by FONLL 
(pt^ = 4.5 GeV/c) and the highest possible crossing point allowed by the calculation (p^^ = 
10.5 GeV/c). The latter results in the weakest contribution of electrons from b quark decay 
to the total non-photonic electron yield. As shown in Fig. 5, whatever the assumed crossing 
point, the effect of the Ac/ D enhancement remains visible. It leads to a decrease of the 
non-photonic electron Raa of about 10(25)% for a crossing point at p^^ = 4.5(10.5) GeV/c. 

In addition to the A^/D enhancement addressed in this work, one could expect an en- 
hancement of the Dg/D ratio due to the strangeness enhancement in heavy ion collisions. 
According to Tab. I, the Df semi-leptonic decay branching ratio is similar to that of 
which represents the main source of non-photonic electrons. An enhancement of Df mesons 
would therefore not affect significantly the Raa of non-photonic electrons. However, as 
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the uncertainty on the measured semi-leptonic decay branching ratio of Df mesons is large 
(Tab. I), the contribution of Df mesons to the non-photonic electron yield and consequently 
to the non-photonic electron Raa cannot be estimated precisely. From the theoretical side, 
according to the Spectator Model for charm mesons decay, the semi-leptonic decay widths for 
the different charm mesons should be equivalent [24]. Knowing charm meson lifetimes [18] 
and branching ratios for D° and D"^ (see Tab. I), one can estimate the Df branching ratio 
to 8.2 ± 0.2% which appears to be consistent with the measured value. 

In summary, we have shown that an enhancement of the Ac/D ratio in heavy ion collisions 
has important consequences on the nuclear modification factor of non-photonic electrons. 
Such an enhancement, which has recently been predicted by the coalescence model and which 
has already been measured for non-strange and strange hadrons, would significantly reduce 
the Raa of non-photonic electrons at intermediate Pf This is a consequence of the smaller 
semi-leptonic decay branching ratio of charm baryons compared to that of charm mesons and 
of the softer decay lepton spectrum from charm baryons compared to that of charm mesons. 
In the most realistic situation investigated in the present work the enhancement leads to an 
additional non-photonic electron suppression of 10 — 25% (with respect to the suppression 
observed without charm baryon/meson enhancement). This suppression can even be larger 
in case of a weaker bottom contribution to the non-photonic electron spectrum. We conclude 
that it is therefore premature to interpret the non-photonic electron Raa data before a 
possible enhancement of the A^/ D ratio is measured experimentally. Heavy quark energy 
loss can be studied in a much cleaner way via the nuclear modification factor of exclusively 
reconstructed charm hadrons. Such measurements should be possible with the RHIC-II 
experiments [25] and with the ALICE experiment at the LHC [26]. We finally note that the 
Ac/ D enhancement can possibly infiuence the elliptic fiow of non-photonic electrons as well. 
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TABLES 

TABLE I. Inclusive decay branching ratio (BR) of charm hadrons into e + anything [18] and 
yield (A^) of charm hadrons (in 47r) in pp collisions at ^/s = 200 GeV from the present PYTHIA 
simulation using input parameters as described in [21]. The total cross-section for charm production 
is normalized to the experimental value obtained in [21]. N^^ and A?^^ include primarily produced 
Ac and Ac as well as those from Sc and Sc decay. ei\f^ {BR) is the contribution to the uncertainty 
of the total electron yield due to the uncertainty on the particle BR. 



Hadron 


D+ D- 


DO 


Dt 


D7 


Ac Ac 


BR (%) 


17.2 ±1.9 


6.71 ±0.29 




811 


4.5 ± 1.7 


(xlO-^) 


3.00 3.07 


9.31 9.85 


1.82 


1.60 


1.23 0.85 


eM^BR) (%) 


1.08 1.10 


1.31 1.39 


4.41 


3.58 


1.51 1.04 
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TABLE II. Main differences between the approach presented in [17] and this work. See text 
for more details. 





i-^ ' J 


this work 


Ac/ D versus pt in Au+ Au coUisions 


as A/K° data 


Gaussian 


Ac/ D versus pt in pp collisions 


as A/X^ data 


PYTHIA 


Maximum of the Ac/i^ enhancement 


~ 1.7 at pt ~ 3 GeV/c 


~ 0.9 at pt ~ 5 GeV/c 


Energy loss 


hadron shape scaling 


[23] 


Electrons from bottom decay 


no 


yes 
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FIGURES 
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FIG. 1. Invariant differential cross-section of non-photonic electrons (dots) measured in pp 
collisions at ^/s = 200 GeV [13]. The dashed curves show the prediction from FONLL calcula- 
tions [22]. The solid curve shows the result of the PYTHIA simulation as described in the text. 
The simulated spectrum is normalized from the integration of the measured spectrum in the range 
1.4 <pt<4. GeV/c. 
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FIG. 2. Transverse momentum dependence of the charm baryon/meson ratio (squares) and 
decay electrons from charm baryons over decay electrons from charm mesons (dots). The results 
are obtained from the P YTHIA simulation described in the text for pp collisions at s/s = 200 GeV. 
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Our study, enhancement factor C=12 
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FIG. 3. Nuclear modification factor of non-photonic electrons (dots) measured in central 
(0 — 10%) Au+Au collisions at ^snn = 200 GeV [5]. The crosses correspond to the results 
of the simulation described in the text for a Ac/D enhancement factor of 12. 
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PHENIX data 
* c, rad energy loss 
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FIG. 4. Nuclear modification factor of non-photonic electrons (dots) measured in central 
(0 — 10%) Au+Au collisions at ^/sJ/n = 200 GeV [5]. The symbols show the result of the simulation 
described in the text. 
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• PHENIX data 

A c+b, rad+coll loss, P^''=4.5 GeV/c 
X c+b, rad+coll loss, P^''=4.5 GeV/c, C=12 J 
c+b, rad+coll loss, P^''=10.5 GeV/c 
c+b, rad+coll loss, P^''=10.5 GeV/c, C=12 
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FIG. 5. Nuclear modification factor of non-photonic electrons (dots) measured in central 
(0 — 10%) Au+Au collisions at ^/sJ/n = 200 GeV [5]. The symbols show the result of the simulation 
described in the text. 
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